arXiv: 1502.00714vl [cs.IT] 3 Feb 2015 


Exploiting the Preferred Domain of FDD Massive 
MIMO Systems with Uniform Planar Arrays 


Junil Chod, Taeyoung Kind, David J. Loved and Ji-yun Seod 

^School of Electrical and Computer Engineering, Purdue University, West Lafayette, IN 
Email: choi215@purdue.edu, djlove@purdue.edu 

^Communications Research Team (CRT), DMC R&D Center, Samsung Electronics Co., Ltd., Suwon, Korea 
Email: ty33.kim@samsung.com, jiyun.seol@samsung.com 


Abstract —Massive multiple-input multiple-output (MIMO) 
systems hold the potential to be an enabling technology for 5G 
cellular. Uniform planar array (UPA) antenna structures are a 
focus of much commercial discussion because of their ability to 
enable a large number of antennas in a relatively small area. 
With UPA antenna structures, the base station can control the 
beam direction in both the horizontal and vertical domains 
simultaneously. However, channel conditions may dictate that 
one dimension requires higher channel state information (CSI) 
accuracy than the other. We propose the use of an additional 
one bit of feedback information sent from the user to the base 
station to indicate the preferred domain on top of the feedback 
overhead of CSI quantization in frequency division duplexing 
(FDD) massive MIMO systems. Combined with variable-rate 
CSI quantization schemes, the numerical studies show that the 
additional one bit of feedback can increase the quality of CSI 
significantly for UPA antenna structures. 

I. Introduction 

Future wireless cellular systems are expected to deploy 
a large number of antennas, e.g., IOs-IOOs of antennas, at 
the base station □-10. This trend of deploying a large 
number of antennas is now well known as massive multiple- 
input multiple-output (MIMO) systems. After pioneering work 
done by Marzetta in £ 9 , many follow up works have been 
dedicated to verifying the benefits, potential uses, drawbacks, 
and limitations of massive MIMO systems. We refer to □, 
0 and references therein for details. 

Most of the previous works on massive MIMO focused 
on time division duplexing (TDD) to circumvent downlink 
channel state information (CSI) estimation and quantization 
problems. Frequency division duplexing (FDD) is extremely 
challenging to implement with massive MIMO because most 
previous solutions for CSI estimation and feedback design 
become impractical as the number of antennas grows large. 
However, most current wireless cellular systems are based 
on FDD, and backward compatibility is crucial for advanced 
wireless communication technologies. Thus, it is of great 
interest to solve CSI estimation and quantization problems for 
massive MIMO systems. 

There are some existing works dedicated to making FDD 
massive MIMO practical. Low complexity CSI quantization 
schemes are developed in m-m where GDI specifically 
focused on backward compatibility with the 3GPP LTE- 
Advanced standard. The schemes in 0, 0, iTTOl are variable- 


rate CSI quantization techniques that adaptively control the 
feedback overhead of CSI quantization. For the channel 
sounding problem, tiD- m showed that training overhead 
can be significantly reduced by adapting the training signals 
using knowledge of the long-term channel statistics. These 
past works, however, have not considered practical antenna 
structures of massive MIMO systems with a few exceptions 
in GO, £D that consider uniform planar arrays (UPAs), also 
sometimes referred to as uniform rectangular arrays (URAs). 

Interest in UPAs for massive MIMO deployment is growing, 
mainly because of the ability of a UPA to house a large number 
of antennas in a small area. Massive MIMO with a UPA is 
sometimes referred to as three-dimensional (3D) MIMO 1 l~3ll— 
IE9 because the base station can control the beam direction 
in both the horizontal and vertical domains simultaneously 
using the two-dimensional structure of a UPA. In this case, the 
channel might require more CSI accuracy for the horizontal or 
vertical domain depending on the scenario. Thus, we propose 
to use one additional bit of feedback to achieve the full benefit 
of variable-rate CSI quantization schemes in UPA scenarios. 

The additional one bit of feedback explicitly indicates 
the preferred domain between the horizontal and vertical 
domains. It is important to point out that the concept of 
the preferred domain in CSI quantization is novel because 
current wireless communication standards such as 3GPP LTE 
and LTE-Advanced only focus on one-dimensional antenna 
array structures tm With one-dimensional antenna array 
structures such as a uniform linear array or dual-polarized 
linear array, the user only sees the horizontal domain beam 
pattern. The key idea of the proposed scheme is that the user 
re-indexes the channel elements to be quantized depending on 
the domain that needs more precise quantization. Simulation 
results show that the additional one bit of feedback with 
appropriate CSI quantization schemes can increase the quality 
of CSI significantly for massive MIMO with UPA antenna 
structures. 

The paper is organized as follows. We explain our system 
model and discuss CSI quantization techniques for UPA struc¬ 
tures in Section HI] In Section Hill we propose our one bit of 
additional feedback idea. Numerical studies that evaluate the 
proposed idea are shown in Section llVl and conclusions follow 
in Section [V] 


II. System Model and CSI Quantization 

We describe our system model first. Then, we explain 
problems and possible solutions of CSI quantization for UPA 
antenna structures. 


A. System Model 

We consider a multiuser (MU) MIMO block-fading channel 
where the base station with N p = N v x Nh antennas serves 
I\ single antenna users simultaneously. The base station is 
equipped with N v rows and Nh columns of UPA antennas 
with antenna spacing d\ for the vertical domain and r/ 2 for 
the horizontal domain. We define the (fc,Z)-th antenna as the 
antenna element in the /c-th row and the Z-th column of a 
UPA, which results in the ( k , Z)-th antenna corresponding to 
the l + Nh{k— l)-th element of h[m]. Fig.Q]shows an example 
of 4 x 8 UPA structure with antenna indexing. 

Assuming equal power allocation, the received signal of the 
user i at the m-th fading block is written as 


Vi[m] = W -^hf [m]wi[m]si[m] 

V A 


K 

- y 

I< ^ 

u=l,u^i 


hf [m]w u [to] s u [to] + [to] 


where hj[m] G C Np , w,[to] G C Np , and Zi[m] ~ CAf(0, 1) 

are the channel vector, the unit norm beamforming vector, 
and the complex additive white Gaussian noise of the z-th 
user, respectively. ,s, [m] is the data signal of the z-th user with 
E [sj [&]] = 0 and 7? [|sj[m]| 2 ] = 1, and p is the transmit 
signal-to-noise ratio (SNR). 

Note that w, [to] is a function of the quantized CSI received 
from all users, i.e., hj[m] for i = 1,..., K. We rely on well- 
known zeroforcing beamforming (ZFBF) for w,; [rn] H 17H . Let 
H[to] be the composite matrix of h, as 


H[to] 


hi [m] ■ ■ ■ Yik [to] 


( 1 ) 


Then, the composite ZFBF precoding matrix W[to] is given 
as _i 

W[m] = H[m] ^HfTOj^Hfrn]^ , 
and the beamforming vector for user i becomes 


Wj [to] = W (;ii )[m] 


where A r. j.) is the fc-th column of the matrix A. 

The instantaneous signal-to-interference-noise ratio (SINR) 
of the z-th user is given as 


SINR, [to] 


| [m]wj [to] 
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| [to] w u [to] 
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and the corresponding sum-rate of the TO-th fading block is 
written as 


K 

Rsum[m] = y log 2 (1 + SINR„[m]) 

If—1 
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Fig. 1: An example of UPA antenna structure of N v = 4 and 
Nh = 8 with antenna element indexing. 



assuming Gaussian signaling. 

We consider spatially and temporally correlated channels in 
this paper. Thus, we model the channel vector of the i-th user 
as 

hi [to] = iphi [to - 1] + yj 1 - rjfRf g, [to] (2) 


where 0 < r/, < 1 is the temporal correlation coefficient, 
Ri = E [hj [to] hj [to] h ] is a spatial correlation matrix, and 
g, [to] ~ CAf(0, 1) is an innovation process with i.i.d. Rayleigh 
fading components. At m = 0, we have h, [0] = R,? g,; [0]. 

For the spatial correlation matrix R,, we adopt the UPA 
spatial correlation model from 0, where the correlation 
between the (fc, Z)-th and (p, q)-th antenna elements is given 
as 


[R-i](fe,I),(p,g) 

with the variables 


nl + (sin i>^)CT^) 2 ■ 


(3) 


D\ = ( p _fe) C os6l ie -i(f i 2^i) 2 (p_fc) 2 s i n 2 0^ 

D l 2 = sin6»i, 

D 3 = l)cos0i, 

D\ = i dld ?(p - k )(q - 0 sin(20i) , 

D l = ( D l) 2 ((sin^i)crQ 2 + 1, 

D 6 = -D 4 (( sin< Z !, i) cr i ) 2 +COS^i, 

D\ = ( D 3 i ) 2 cos 2 <j)i - (I? 4 t) 2 ((sin fa)^) 2 - 2D\ cos </>i 


where A is the carrier frequency wavelength, fa is the mean 
horizontal angle-of-departure (AoD), (/, is the mean vertical 
AoD, <Ji is the standard deviation of horizontal AoD, and fa 
is the standard deviation of vertical AoD of the z-th user. The 
boresight angle of fa and 9i are both n/2. 


B. CSI Quantization for UPA 

Because CSI quantization is performed independently at 
each user, we drop the user index i in this subsection. To 
obtain a reasonable CSI quality, we assume that the total CSI 
quantization overhead Rcsi is given as 


Rcsi = BN P 
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( 4 ) 


where B is the quantization bits per antenna element. If we 
rely on the conventional approach of using A' r ,-dimensional 
vector quantized codebooks for CSI quantization, the computa¬ 
tional complexity that grows exponentially with f?csi becomes 
a serious problem for large N p 0, which is typically the case 
of UPA antenna structures. 

There are several ways to solve the CSI quantization com¬ 
plexity issue for UPA antenna structures. It was shown in (§1 
that the spatial correlation matrix R can be approximated a.^j 


R « R , 0 R h 


(3) 


where 0 is the Kronecker product and 

[Ru]fc,p = e -i(f^) 2 (p-fe) 2 sin 2 0 e j2^ ( P -fe )co S 0 ; 

f £?3 cos 2 (sin <ft)cr) 2 .£2 COS 4 > 

[R/i]z, 9 = — f =e 3355 e J n5 

V -^5 

are the spatial correlation matrices of the vertical and horizon¬ 
tal domains, respectively. 

The approximation in (0 suggests to quantize the CSI of 
horizontal and vertical domains separately. Let H[m] be the 
matrix consists of the elements of h[m] all 


^Nh + l-^Nh] M 


h hV4tVh —l) + l:JV„iV h ] M 


where a[ rai:ri2 ] is the sub-vector of a consists of the ?7i-th to 
7(2-th entries. Then the user can quantize the CSI of horizontal 
and vertical domains as 


c/i.optM = argmax 11 H[to]c*. 11 , 
c h&C h 

c„, op tM = argmax ||c^H[m]|| 
c„ec„ 

with Bh- bits, ./V^-dimensional codebook C/, and /I,,-bits, N v - 
dimensional codebook C v with a constraint Bh + B v = i?csi 
(or Rcsi + 1 taking the proposed one bit of additional 
feedback into account). Once the user feeds back c^ opt [m] 
and 0 pt [to], the base station can reconstruct the quantized 
CSI as 

hH = c„ i0pt [m] 0 c* h opt [m}. 


1 In (SI- the approximation is given as R. m 10, (8> R , because 0 indexed 
antenna elements vertically. 

2 The /)-th element of H[m] corresponds to the (k. /)-th antenna 
element of the UPA antenna structure. 


We dub this approach the Kronecker-product approach and 
consider it as the baseline of performance evaluation^ 

On the other hand, we can quantize 1 i[to] in a block-wise 
manner. For example, we can quantize L <C N p dimensional 
sub-channel vectors separately using BL bits for each sub¬ 
channel vectorQ The problem is that the codeword selection 
criterion of using an N p -dimensional codebook cannot be 
decomposed into N p /L sub-problems because Q 

argmax | h H [m] c | 2 
c ec 

N p /L 

argmax Ei hf L(n-l)+l:In] H C [L(n-l)+l:Iii] P 

c6C n—1 

assuming L divides N p . However, we can transform the 
problem as in 0. Thus, we can quantize each separate block 
of 1 i[to 1 with a small number of BL bits in a noncoherent 
mannerjj which can reduce the complexity significantly. 

Taking this transformed problem into account, the trellis- 
extended codebook (TEC) and trellis-extended successive 
phase adjustments (TE-SPA) that quantize h[m] in a block- 
wise manner have been proposed in Col. To be specific, 
TEC is a base codebook that quantizes h[0], and TE-SPA is 
a differential codebook m-m that exploits the temporal 
correlation of channels and quantizes h[m] for to > 1 using 
the previously quantized CSI h[m — 1]. Another advantage 
of TEC and TE-SPA is that they support variable-rate CSI 
quantization, which makes it easy to adapt feedback overhead 
depending on requirements. In this paper, we assume the one 
bit of additional feedback idea that will be explained in the 
next section is combined with TEC and TE-SPA. However, the 
proposed one bit of additional feedback idea can be applied 
to any kind of block-wise CSI quantization schemes. 

III. One Bit of Additional Feedback 

Note that a UPA is a two-dimensional antenna structure 
and can control the beam direction not only of the vertical 
domain but also of the horizontal domain. Moreover, the 
approximation in 0 shows that the CSI of one domain might 
need more accurate quantization than the other. Thus, we 
propose to have one bit of additional feedback, which indicates 
the preferred domain, from the user to the base station. 

'Although not exactly the same, similar CSI quantization techniques as the 
Kronecker-product approach have been proposed in 0, 0. 

4 Note that L is the system parameter that needs to be optimized to minimize 
CSI quantization error with fixed Bcsi> which is out of scope of this paper. 

5 Note that i/i is only an auxiliary variable during the optimization process, 
and the user does not need to feed back the information of -i/i to the base 
station. We refer to ED for details. 








We drop the user index i and the fading block index m to 
simplify notation. We assume N v and Nh are multiples of L. 
We further assume that h is quantized in a block-wise manner 
with the block size L. Using the L-dimensional codebook T , 
we first generate a candidate quantized CSI hi as 

i T 


hi = 


fl . . . fl 

r l,l T l. 


N p /L 
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(a) The channel vector h is quantized with hi. 


(fl,l) ■ ■ ■ > ft ,n p /l) — 

N p /L 

argmin min V* ||h[ L ( n _ 1)+1 . Lrl | - e^fi^ll . 
fi.,£7 i>£[0,2n) ^ 

This problem can be solved by TEC, TE-SP/fJ or any other 
block-wise CSI quantization schemes. We assume the elements 
of T are properly normalized to have ||hi|| 2 = 1. 

In addition to hi in ©, we generate another candidate 
quantized CSI I 12 as 


h 2 



f T 

L 2,N p /L 


(5) 


with 

(f 2 ,t, ■ • •)f 2 ,JV p /Z,) = 

N p /L 

argmin min Y] hr L („_ 1)+1 . Ln ] - e H ’t 2 ,n 

Ve[0, 2 7T) II 1 ^ ' 


where h is the rearranged channel vector of h. To explicitly 
describe h, we define two functions 


c(n) = 


nL 


m(n) = 



1) mod 



With these functions, h is given as 


h = 


hf ..-h T 


N p /L 


T 


where 


hn — \p j c{ri)+Lm{n)Nh^c{n)+{Lra{n)+l)Nh'> 

i ^c(n)+(L(m(n)+l)-l)Nh\ 


for n = 1,..., N p /L. The final quantized CSI h is then given 
as 



if |h H hi| 2 > |h H h 2 | 2 
if |h H hi| 2 < |h ff h 2 | 2 


Note that we can adopt an arbitrary L-dimensional code¬ 
book T when we solve (Q} and ©, and it is well known 
that standard codebooks such as DFT or LTE codebooks have 
good beam directivity when antenna spacing is small. Thus, it 
would be a good choice to adopt DFT or LTE codebooks as 
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(b) The rearranged channel vector h is quantized with I 12 . 

Fig. 2: The conceptual explanation of the proposed idea with 
4x8 UPA antenna structure and L = 4. The red box is 
quantized with fi >n in (a) and f 2j „ in ( b ) for n = 1,..., N p /L 
where the numbers in the red boxes correspond to n. 


T for block-wise CSI quantization schemes for UPA antenna 
structures. 

Now, we conceptually explain the proposed idea usin<£] 
Fig. [2] It might be beneficial for the user to quantize one 
domain more precisely than the other depending on the sce¬ 
nario, e.g., user location or user movement. If the user needs 
to quantize the horizontal domain more precisely, than the 
user should quantize the channel horizontally using structured 
codewords. On the other hand, the user needs to quantize 
the channel vertically with structured codewords if the user 
needs more accurate CSI for the vertical domain. With these 
observations, the candidate CSI hi quantizes the horizontal 
domain of the channel more precisely as in (a) of Fig. [2] 
while h 2 quantizes the vertical domain of the channel more 
accurately as in (b) of Fig. [2] By comparing hi and h 2 , the 
user can determine the more important domain for the user and 
select more accurately quantized CSI. The additional one bit 
of feedback indicates which candidate CSI is selected between 
hi and h 2 . 

Remark 1: If we rely on the conventional approach of 
using N p -dimensional vector quantized codebook for CSI 
quantization, then the proposed one bit of additional feedback 
has eventually the same effect of using a f?csi + 1-bits vector 
quantized codebook. However, block-wise CSI quantization 
schemes are not able to exploit additional few bits (one bit 
in our propopsed idea) of feedback in this manner. Thus, 
the proposed idea is particularly suitable to block-wise CSI 
quantization schemes. 


6 Using a trellis structure, we can have a codebook expansion effect in TEC 
and TE-SPA, e.g., 8 codewords are available for T while only allocating 2 
bits per sub-problem in CD- We refer to GSD for more details. 


7 A figure similar to Fig. [2 was shown in ED; however, the figure in (25] 
was to propose downlink training design instead of CSI feedback. The concept 
of the preferred domain for CSI feedback in this paper is new. 



























Fig. 3: Sum-rate (bps/Hz) according to N p = N v x Ay, anten¬ 
nas. The feedback overhead of TEC from [ 101 is Bqsi = \N p 
bits while that of other two schemes is Bq si + 1 bits. 

Remark 2: The channel vector rearrangement should be de¬ 
pendent on antenna indexing and block-wise CS1 quantization 
schemes. However, the proposed idea can be applied to any 
scenario without difficulty. 

IV. Simulation Results 

We evaluate the proposed idea with Monte-Carlo simu¬ 
lations in this section. We assume the base station serves 
K = 10 users simultaneously with ZFBF prccodingQ We 
adopt the channel model in (0 with the spatial correlation 
matrix in <0. Considering a practical UPA structure and cell 
layout, we set the vertical antennna spacing to <1\ = 0.9A, 
horizontal antenna spacing to d ,2 = 0.5A, randomly uniformly 
generate the mean horizontal AoD fa in [7t/6, 57t/6), the mean 
vertical AoD Oj in [tt/ 12, tt/ 8), and the standard deviations of 
horizontal AoD ay and vertical AoD both in [7t/18, 7t/12) 
independently for all users in each channel realization. The 
transmit SNR p is set to HMB for all simulations. 

We first compare the MU-MIMO sum-rate of three schemes, 
i.e., the Rronecker-product approach explained in Section 
III-BI TEC from fTO) , and TEC with the proposed one bit 
of additional feedback. We set L = 4 and adopt the L- 
dimensional DFT codebook that satisfies B = 1/2 bit per an¬ 
tenna quantization for TEC. The Rronecker-product approach 
is based on two N v and ^-dimensional DFT codebooks of 
sizes B v = jN p and B^ = jN p + 1 bits. The feedback 
overhead of TEC is given as Be si = \N p while that of the 
other two schemes is Bcsi + 1- We only consider the m = 0- 
th fading block in this simulation because all three schemes 
do not exploit any temporal correlation, . 

In Fig. [3] we plot the sum-rates of the three schemes 
according to different combinations of N p = N v x Ay, UPA 

8 To perform ZFBF precoding properly, we only consider the channel 
realizations when the composite quantized channel matrix H[m] given in 
0 is full rank. This can be thought as a very simple scheduling scheme. 



Fig. 4: Sum-rate (bps/Hz) according to the fading block index 
to with 8x8 UPA antennas. TEC quantizes CSI at to = 0 and 
TE-SPA quantizes CSI at to > 1. The feedback overheads of 
all three schemes are the same as in Fig. [3] for all to. 


antennas. It is clear that the proposed one bit of additional 
feedback gives around 0.5 to 0.8 bps/Hz gain of sum-rate 
depending on antenna configurations. Moreover, TEC with 
the proposed one bit of additional feedback is comparable to 
the Rronecker-product approach in the 4x8 UPA antenna 
scenario and outperforms the Rronecker-product approach in 
the 6x8 and 8x8 UPA antenna scenarios. Although we 
could not simulate higher dimensions of UPA because of the 
complexity issue of the Rronecker-product approach (e.g., a 
20 bit codebook is needed for the Rronecker-product approach 
when N p = 10 x 8), we expect that TEC with the proposed 
one bit of additional feedback would keep outperforming 
the Rronecker-product approach in larger UPA dimensions as 
well. 

If Fig. [4] we plot the sum-rates of three schemes according 
to the fading block index to with 8x8 UPA antennas. We adopt 
Jakes’ model to generate the temporal correlation coefficient 
as l[26l 

V = J 0 (27 Tf D r) 

where Jo(-) is the zeroth-order Bessel function, //; is the 
Doppler frequency, and r is the channel instantaneous interval. 
With 2.5GHz carrier frequency, 3km/h user velocity, and 5ms 
channel instantiation interval, we have p = 0.9881. Note 
that the Rronecker-product approach does not exploit temporal 
correlation of channels while other two schemes exploit the 
temporal correlation by quantizing CSI with TE-SPA flOl for 

TO > 1. 

The figures clearly show that the two block-wise CSI quan¬ 
tization schemes outperform the Rronecker-product approach 
as m increases. We might be able to take advantage of the 
temporal correlation even for the Rronecker-product approach; 
however, to the best of our knowledge, no such scheme 
has been developed yet. Note that the proposed one bit of 



















additional feedback significantly increases CSI quality for all 
to, which verifies its effectiveness. 

V. Conclusion 

We proposed to have one bit of additional feedback on top of 
CSI quantization overhead from the user to the base station for 
FDD massive MIMO systems in this paper. The proposed one 
bit of additional feedback indicates the preferred domain of the 
user’s channel, which is a new concept in CSI quantization. 
This concept has become available because of the use of 
UPA antenna structures that are getting more interest due to 
massive MIMO systems. The proposed idea is well suitable 
to block-wise CSI quantization schemes which are shown to 
be appropriate for UPA antenna structures. The simulation 
results showed that the proposed one bit of additional feedback 
idea can significantly increase the quality of quantized CSI in 
practical UPA channel conditions. 
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